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Q , Abstract 

\D 

fS| ' We present a one-loop calculation of the static potential in the SU(2)-Higgs model. The 

^O , connection to the coupling constant definition used in lattice simulations is clarified. The 

consequences in comparing lattice simulations and perturbative results for finite temperature 

0^ ■ applications are explored. 
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D , 

^. : 1 Introduction 

>• 

^ - The observed baryon asymmetry of the universe was eventually determined at the electroweak phase 
H ,. transition |jl[]. The most straightforward method to study this phase transition is to use resummed 
perturbation theory (cf. e.g. 0, |^, Q). In the low temperature Higgs phase the perturbative approach 
is expected to work well, however, it is not able to describe the high temperature symmetric phase, 
which has serious infrared problems in perturbation theory. Since the determination of thermodynam- 
ical quantities at the critical temperatures is based on the properties of both phases, non-perturbative 
techniques are necessary for a quantitative understanding of the phase transition. 

One very succesful possibility is construct an effective 3-dimensional theory by using dimensional 
reduction, which is a perturbative step. The non-perturbative study is carried out in this effective 
3-dimensional model (see e.g. |^ and references therein). Analytical estimates are confirmed by 
numerical results and relative errors are believed to be at the percent level. 

Another approach is to use 4-dimensional simulations. The complete lattice analysis of the 
Standard Model is not feasible due to the presence of chiral fermions, however, the infrared problems 
are connected only with the bosonic sector. These are the reasons why the problem is usually studied 
by simulating the SU(2)-Higgs model on 4-dimensional lattices, and perturbative steps are used to 
include the U(l) gauge group and the fermions. Finite temperature simulations are carried out 
on lattices with volumes Lt • L^g, where Lt <C Lg are the temporal and spatial extensions of the 
lattice, respectively. The lattice spacing is basically fixed by the number of the lattice points in the 



temporal direction {Tc = l/{Lta), where Tc is the critical temperature in physical units); therefore 
huge lattices are needed to study the soft modes. This problem is particularly severe for Higgs boson 
masses around the W mass, for which the phase transition is weak and typical correlation lengths are 
much larger than the lattice spacing. In this case asymmetric lattice spacings are used, in particular 
the spatial lattice unit is approximately four times larger than the temporal one . 

Despite the fact that the two approaches (perturbative and lattice) are systematic and well- 
defined, it is not easy to compare their predictions. The reason for this is that in lattice simulations 
the gauge coupling constant is determined from the static potential, whereas in perturbation theory 
the MS scheme is used. The main goal of this paper is to perform as perfect a comparison as possible, 
by determining the MS gauge coupling constants, which correspond to the different lattice results.^ 

The paper is organized as follows. Section 2 contains the one-loop static potential of the SU(2)- 
Higgs model using the MS scheme in the Feynman gauge. Section 3 relates the continuum version 
of the lattice gauge coupling constant to the MS coupling. In Section 4 the detailed comparison of 
lattice and perturbative predictions are presented. Section 5 summarizes our results. 

2 Calculation of the one-loop static potential 
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The one-loop static potential was calculated long ago in quantum chromodynamics p|, 
even the full two-loop result was published recently |Tl[|. The calculation is based on the same 
principles and techniques in the case of the SU(2)-IIiggs model. One calculates rectangular Wilson 
loops of size r xt. The logarithm divided by —t gives the potential at distance r in the t —^ oo limit. 
Our calculation was performed in the MS scheme and the Feynman gauge but the result is gauge 
independent, as it should be for a physical observable. The relevant graphs are shown in Fig. 1. 
Other graphs, giving vanishing contributions in the Feynman gauge and are not shown in Fig. 1. 
Solid lines represent the heavy quark (antiquark) propagator, while wavy lines the vector boson 
propagator. External heavy quark (antiquark) propagators are not shown in the figure. The one- 
loop corrected vector boson propagator contains scalar and ghost contributions as well. The result 
can be conveniently given in momentum space. One obtains 
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^ During the write-up of our results, prior to us, a similar, independent ealculation for the gauge coupling constant 
was presented by M. Laine H, who compared 4-dimensional and 3-dimensional results, too. Using his convention for 
the renormalized gauge coupling, which is a special case of our definition, the two results agree. However, as it will 
be discussed later, our definition for the perturbative gauge couplings is closer to the actual lattice definitions. 
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Figure 1: Graphs giving nonvanishing contributions to the static potential 
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where k"^ denotes the square of the three-momentum k, Mh the Higgs mass and Rhw = Mh/Mw- 
The function F is defined as 
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As it can be seen, our result does depend on the renormalization scale fi and it fully agrees with that 
of M. Laine 0. 

Eq. (P has to be Fourier transformed into coordinate space. We applied the brute force method 
performing numerical integration. As a check, we compared our results with various pieces of the 
partly analytic calculation in |^ for the derivative of the potential (with respect to distance). The 
agreement is excellent. 

Our result is presented in Figs. 2 and 3, where the various parts of the one-loop correction to the 
potential are plotted. We define 
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where M^r = My/ — 6Mw, with 6Mw the one-loop mass correction. Since 6Mw is scale dependent, 
so is M{y. A and B are functions of the distance r and Rhw = Mh/Mw- We choose Mw = 80GeV. 
Fig. 2 shows the dependence of A and B on the dimensionless distance rMw for Rhw = 0.8314 
(corresponding to the end point of the first order finite temperature phase transition |]12|); while 
Fig. 3 shows the Rhw dependence for r = M^^. 



3 Relation of the continuum version of the lattice coupling 
constant definition to the MS coupling constant 

Since we wish to compare results of lattice simulations and continuum perturbation theory calcu- 
lations, it is an essential point to define the SU(2) gauge coupling in the same way in both cases. 
However, in continuum perturbation theory the MS running coupling constant at a given renormal- 
ization scale is more natural (as used in Eqs. (|l],0), too), while in lattice simulations other definitions 
are applied. Therefore we have to establish the relation between the coupling constants. 
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Figure 2: The coefficients of (/^/(IGtt^) — curve A — and of 5''*/(167r^) log(yU^/M^) — curve B — of the 
static potencial Eq. (|) as a function of distance times W mass. i?jyvF=0.8314. 
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Figure 3: The coefficients of g'^ j {\&t\:'^^ — curve A — and of 5''^/(167r^) log(/i^/M^) — curve B — of the 
static potencial Eq. (|) as a function of Rhw- The distance is M^^. 



The lattice definition of tlie coupling constant is given in [^ . Note that we are using the local 
version below. For the reader's convenience we recall this definition (inspired by |I^). 
First rectangular Wilson loops of size (r,t) are measured. Extrapolating to large t and dividing 
the logarithm by —t one gets the static potential in the t ^ oo limit as a function of r. The 
nonperturbative lattice static potential is fitted by a finite lattice version of the Yukawa potential 
with four parameters (for details cf. |T^). One of these parameters is the mass in the exponential of 
the Yukawa potential, which is usually called the screening mass. The gauge coupling at distance r 
is defined as the ratio of the discrete r derivative of the lattice simulated nonperturbative potential 
and the discrete derivative of the tree-level lattice Yukawa potential normalized by the square of 
the tree-level coupling and with the mass parameter Miattice identified with the screening mass. In 
practice 5'iattice(Mattice) is determined and is called the local renormalized gauge coupling constant 
on the lattice. The lattice results at various Higgs masses are collected in Table 1. Data are from 
§,0,0, and [0. 

To follow the above procedure in the case of the continuum perturbative determination of the 
renormalized gauge coupling, we performed a fit of the one-loop potential with a tree-level Yukawa 
potential plus a constant term. The parameters of the fit are the coupling constant, the mass in the 
exponent (perturbative "screening mass" Mgcreen) and the constant. For the various values of Higgs 
mass we performed the fit in the same r range as used in the lattice studies and took the errors of the 
fitted function to be proportional to the errors of the potential obtained in lattice simulations. g]i{r) 
is then determined by taking the ratio of the derivatives with respect to r of the one-loop potential 
and the tree-level potential normalized by the square of the tree-level coupling, i.e. we have 
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with V{r) given by Eq. (H), Cp = 3/4, and Mscreen 

for distances satisfying Mgcreen — V''" = ^(fi'^) ^^ can put Eq. 



obtained from the fit. Since Mg, 
into the form 
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C and D are functions of Rhw and Mgcreen, 
80GeV. 



their values are tabulated in Table 2 for Mo, 



M, 



w 



Rhw 


.2049 


.4220 


.595 


.8314 


Tc (GeV) 


38.3 


72.6 


100.0 


128.4 


Miattice (GeV) 


84.3(12) 


78.6(2) 


80.0(4) 


76.7(24) 


^littice(M-l) 


.5630(60) 


.5788(16) 


.5782(25) 


.569(4) 


M3,,een (GeV) 


74.97 


80.44 


80.70 


81.77 


9UT.) 


0.540 


0.592 


0.585 


0.570 


cr^iT^) 


0.589 


0.589 


0.579 


0.562 



Table 1: Various quantities calculated for values of Rhw used in lattice simulations. For more 
explanation see the text. As usual the numbers in the parentheses denote the errors in units of the 
last decimals. The errors of the different gauge couplings are dominated by the lattice simulation 
errors (fourth row), therefore we did not indicate them in rows 6 and 7. 



In this procedure we have to choose the gauge couphng in the one-loop potential so that 
gjiiM'^^^^) reproduces the lattice result (third row of Table 1) for the appropriate value of the 
Higgs mass. For our applications (thermodynamical quantities at and around the critical tempera- 
ture Tc of the first order electroweak phase transition) the scale of the one-loop potential is chosen 
to be Tc ~ 2Mh, where Mh is the Higgs boson mass at zero temperature. Thus the gauge coupling 
appearing in the one-loop potential is actually the MS gauge coupling at scale Tc. The MS gauge 
coupling values obtained from this procedure are given in the sixth row of Table 1. 

Another definition of the continuum perturbation theory one-loop "renormalized gauge coupling" 
at distance r is given in |M|. It reads 
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where M is a free mass parameter satisfying M — Mw oc g'^. For this M it is possible to show that 
9R,Laine{^~^) ^au be cxprcsscd in terms of g^^^^Mw), (where Mw is the physical (one-loop) pole 
mass) and all the scale dependence is included in g^j^^Mw)- Assuming M = Mw, the numerical 
difference between this definition and ours is small. However, we believe that it is our definition which 
is the closest conceivable to the local renormalized lattice gauge coupling of [0. In Table 1 (last 
row) we give g^"'^"''^ (Tc) as calculated using Eq. (||), equating gJi^Lainei^w) with the values of the 
lattice simulation results S'latticel lattice) ^^^ using the renormalization group equation to extrapolate 
to the scale T^. 
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-41.54 


-22.19 
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-6.47 


-1.12 


0.5 


-5.66 


1.39 
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-5.23 


2.74 
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-4.98 


3.55 


0.8 


-4.83 


4.06 


0.9 


-4.72 


4.39 


1.0 


-4.65 


4.62 


1.1 


-4.59 


4.78 


1.2 


-4.54 


4.89 


1.3 


-4.50 


4.98 


1.4 


-4.45 


4.98 


1.5 


-4.40 


5.01 



Table 2: The coefficients C and D defined in Eq. (|]) as a function of Rhw- 



4 Comparison of physical observables determined in lattice 
simulations with perturbative predictions 

In the previous section we presented a calculation connecting the renormalized gauge coupling con- 
stant of the MS scheme and g\ obtained from the static potential at different distances. In this 
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section we compare the lattice results and the perturbative predictions for the finite temperature 
electroweak phase transition. Lattice Monte Carlo simulations provide a well-defined and systematic 
approach to study the features of the finite temperature electroweak phase transition. During the 
last years large scale numerical simulations have been carried out in four dimensions in order to clar- 



ify non-perturbative details , |]T2| , |jT3| , |^ . Thermodynamical quantities (e.g. critical temperature, 
jump of the order parameter, interface tension, latent heat) have been determined and extrapolation 
to the continuum limit has been performed in several cases. Nevertheless, it has proven difficult to 
compare the perturbative and the lattice results, because the perturbative approach used the MS 
scheme for the gauge coupling, whereas the lattice determination of the gauge coupling has been 
based on the static potential. The main reason for performing the one-loop calculation of the static 
potential is this kind of comparison. 

In this paper we use the published perturbative two-loop result for the finite temperature ef- 
fective potential of the SU(2)-Higgs model Q]. Note that the numerical evaluation of the one-loop 
temperature integrals gives a result which agrees with the approximation based on high temperature 
expansion within a few percent. The reason for this is that the perturbative expansion up to order 
g"^, X^ corresponds to a high temperature expansion, which is quite precise for the Higgs boson masses 
we studied. It is known that the perturbative loop expansion becomes unreliable for Higgs masses 
above approximately 50 GeV (e.g. resummed perturbation theory fails to predict the end-point of 
the electroweak phase transition, thus it gives a first order phase transition for arbitrarily large Higgs 
boson masses). In the physically relevant range of the parameter space the electroweak phase tran- 
sition can only be understood by means of non-perturbative methods. Therefore it is particularly 
instructive to see quantitatively how the perturbative and the lattice results agree for small Higgs 
boson masses and how they differ for larger ones. 

Since the finite temperature electroweak phase transition is fairly strong for Higgs boson masses 
below 50 GeV, lattices with symmetric lattice spacings were used for Mh ~ 16 GeV, M^ ~ 34 GeV 
and Mfj ~ 48 GeV. The phase transition gets weaker for larger Higgs boson masses, therefore Monte 
Carlo simulations for masses near the W-boson mass are technically difficult. For this parameter 
region different lattice spacings were used in the temporal and the spatial directions. For this type 
of lattice regularization the approach to the continuum limit is somewhat slower; however, even in 
this case it was possible to perform a continuum limit extrapolation for Mh ~ 67 GeV. 

In lattice simulations the gauge coupling constant are determined from the static potential, 
whereas masses are extracted from correlation functions. On the one hand the calculation of the 
previous section connects the gauge coupling definitions between the MS scheme and the scheme 
based on the static potential. On the other hand one can use the zero temperature effective potential 
in order to include the most important mass renormalization effects. The Higgs boson mass obtained 
from the asymptotics of the correlation function corresponds to the physical mass determined by the 
pole of the propagators, i.e. the solution of p^ — M^ = H(p^), where H(p^) is the self-energy. The 
effective potential approach suggested by Arnold and Espinosa approximates H(p^) by H(0) in 
the above dispersion relation. It has been argued that the difference between the two expressions 
is of order g^v"^ {v is the zero-temperature vacuum expectation value), which does not affect our 
discussion. In this scheme the correction to the MS potential reads 



6m^ + --6X + —^\ (7) 





where 

2567r2' 2567r2 V ^ /i 

Here /i is the renormalization scale and My/ is the W-boson mass at T = 0. The above notation 
corresponds to a tree-level potential of the form rri^Lp^ /2 -\- Xip'^/A. Note that this treatment is 




Mh 


16.4(7) 


33.7(10) 


47.6(16) 


66.5(14) 


9l 


0.561(6) 


0.585(9) 


0.585(7) 


0.582(7) 


TJMh 


pert 
nonpert 


2.72(3) 


2.28(1) 


2.15(2) 


1.99(2) 


2.34(5) 


2.15(4) 


2.10(5) 


1.93(7) 


'^+/Tc 


pert 
nonpert 


4.30(23) 


1.58(7) 


0.97(4) 


0.65(2) 


4.53(26) 


1.65(14) 


1.00(6) 





Q/n 


pert 


0.97(7) 


0.22(2) 


0.092(6) 


0.045(2) 


nonpert 


1.57(37) 


0.24(3)* 


0.12(2) 





ojTl 


pert 


0.70(10) 


0.067(6) 


0.022(2) 


0.0096(5) 


nonpert 


0.77(11) 


0.053(5)* 


0.008(2)* 






Table 3: Comparison of the perturbative and the lattice results. The Monte Carlo data are from 
P , |12| , [|1^ , ||15[ (in some cases we have refined the analysis in order to have a more accurate lattice 
prediction). Note that for the mass of the W boson — the dimensionful quantity setting the scale of 
the theory — 80 GeV is used. 



analogous to previous comparisons of the perturbative and lattice results |T^ . 

In m , [|l3l 5 EH 5 EH several observables were determined, including renormalized masses at zero 
temperature [Mr-, ^w'), critical temperatures {Tcj^ jumps of the order parameter (v5+), latent heats 
(Q) and surface tensions (a) for different Higgs boson masses. As usual, the dimensionful quantities 
were normalized by the proper power of the critical temperature. (This convention is adapted in the 
present paper, too.) The simulations were performed on Lt = 2,3,4,5 lattices (Lj is the temporal 
extension of the finite-temperature lattice) and whenever it was possible a systematic continuum 
limit extrapolation was carried out assuming standard 1/a^ corrections for the bosonic theory. 

The statistical errors of these observables are normally determined by comparing statistically 
independent samples. Jackknife and bootstrap techniques were used |T^ and correlated fits were 
performed |jl8[ to obtain reliable estimates of the statistical uncertainties. The systematic errors due 
to finite lattice spacing can be obtained by 1/a^ extrapolation. In cases where it was possible to carry 
out the continuum limit extrapolation we saw that the difference between the Lt = 2 and the Lt = 3 
data was a fairly good estimator of the systematic error. Whenever the data did not make it possible 
to carry out the systematic extrapolation the difference between the Lt = 2 and the Lt = 3 results 
was used to estimate the systematic error. As a conservative estimate we added the statistical and 
systematic errors linearly. For some of the quantities only Lt = 2 data exist. In these cases only the 
statistical errors are listed and an asterisk is used in Table 3 as an indication. A correct comparison 
has to include errors on the parameters used in the perturbative calculation. These uncertainties are 
connected with the fact that neither the Higgs boson mass nor the gauge coupling constant can be 
determined exactly in lattice simulations. Including these errors, the perturbative prediction for an 
observable is rather an interval than one definite value. 

To obtain a better measure of the correspondence between perturbative and nonperturbative 
results, and to incorporate their errors, one introduces "pulls" defined by the expression 



pull 



perturbative mean — nonperturbative mean 
perturbative error + nonperturbative error 



(9) 



The four different pulls at different Higgs boson masses are tabulated in Table 4 and plotted in Fig. 4. 
For the sake of convenience, we used the shorthand Pt = pull of Tc/Mh, P4, = pull of ip+/Tc, 
Pq = pull of Q/T^, and P, = pull of a/T^. 

The quantity which has the smallest pull even for large Higgs boson masses is T^/Mh- A quadratic 
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niH (GeV) 


16.4(7) 


33.7(10) 


47.6(16) 


66.5(14) 


Pt 


4.75 


2.60 


0.71 


0.67 


n 


0.47 


-0.33 


-0.3 


32.5 


Pq 


-1.36 


-0.4 


-1.08 


22.5 


p. 


-0.33 


1.27 


3.5 


19.2 



Table 4: Values of the four different pulls for various Higgs boson masses 
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Figure 4: "Pulls" plotted against the Higgs mass. Arrows indicate values outside the interval [—5, 5]. 



fit was performed to this quantity as a function of Rhw- The result is 



M, 



H 



2.494 - 0M2Rhw + 0.223i?^t^. 



(10) 



In [|12| the end point result for the four-dimensional SU(2)-Higgs model was perturbatively con- 
verted to the Standard Model. In this step the deviation between the two definitions of the gauge 
coupling (f^ljg and the one based on the static potential) was neglected. The estimated uncertainty 
due to this simplification was included into the systematic error of the end point Higgs mass. Using 
the results of the present paper we can also refine the value of the Standard Model end point Higgs 
mass. This is done by perturbatively taking into account fermions and the U(l) factor of the Stan- 
dard Model (for details see the fourth paper of 0). The improvement is established by precisely 
converting the lattice simulation renormalized gauge coupling of the SU(2)-Higgs model into the per- 
turbative g^{Mw)- The new value of the Standard Model end point Higgs mass is 72.1 ± 1.4GeV. 

0] 



This does not deviate much from the old value 72.4 ± 1.7GeV of 



However, the error is smaller, 



since the uncertainty arising from the gauge coupling definitions is eliminated. 

5 Summary 

In this paper we presented the one-loop static potential in the SU(2)-Higgs model. This calculation is 
in agreement with . Using the potential it was possible to connect the gauge coupling constant used 
in finite temperature field theory and in lattice simulations. As expected the numerical difference 
between the two conventions is not that large, it is within a few percent. With this connection 
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we were able to perform a precise comparison between the predictions of perturbative and lattice 
approaches. 

We reanalyzed the existing lattice data and performed a continuum limit extrapolation whenever 
it was possible. The relationship between the two definitions of the gauge coupling constants turned 
out to be marginal, as the lattice data have errors, usually larger then this few percent. The only 
quantity which is measured so precisely that the definition of the gauge coupling constant is essential 
is the ratio of the critical temperature to the Higgs boson mass. As it has been observed already 
for Mh ~ 35 GeV the perturbative value of T^ is larger than in lattice simulations. This sort 
of discrepancy disappears for larger Higgs boson masses. A plausible reason for this fact is the 
convergence of the high temperature expansion used in the perturbative approach. In two-loop 
perturbation theory one uses the high temperature expansion also up to second order, which might be 
inaccurate for the smallest Higgs mass case with temperatures ~ 50 GeV and Higgs field expectation 
values ~ 200 GeV. Nevertheless, the observed differences are on the percent level and they do not 
affect the electroweak phase transition significantly. For small Higgs boson masses (16 and 34 GeV) 
we expect similar differences between lattice and perturbative predictions for other quatities (the 
jump of the order parameter, the latent heat and the surface tension); however, present lattice data 
have too large errors and the differences cannot be seen yet. 

The most dramatic differences appear clearly as we get closer to the end point. The perturbative 
approach gives nonvanishing jump of the order parameter, nonvanishing latent heat and interface 
tension, while the lattice results suggest rapid decrease of these quantities and no phase transition 
beyond the end point. 

Using the results of the present paper we refined the value of the Standard Model end point Higgs 
mass of [12] and obtained 72.1 ± 1.4GeV. 

This work was partially supported by Hungarian Science Foundation Grants under Contract No. 
OTKA-T22929-29803-M28413/FKFP-0128/1997. 
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